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Abstract 

We consider a system of spinless fermions in a nearly one-dimensional cylindrical trap, 
with the Fermi level such that only the lowest-lying states with angular momenta £ = +h, 
0, or — h about the axis of the trap are occupied. We treat the particles in these states as 
comprising three separate Luttinger liquids, with the possibility for forward and backward 
scattering between them. We determine the normal modes in the system in the presence 
of this scattering, finding two sets of linear harmonic fluid modes and one set of modes 
that are self-coupled through a sine-Gordon interaction. 
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The problem of spinless fermionic particles trapped in a cylindrical potential with a 
very high aspect ratio possesses many interesting aspects. If the fermions are constrained 
to move almost exclusively along the direction of the traps' axis, then the system becomes 
effectively one-dimensional in character. This will affect the qualitative nature of the 
system's behavior, since fermions in one dimension have many curious properties. 

In particular, a one-dimensional system of low-temperature fermions exhibits bosoniza- 
tion properties PjJ 12 E| , as the coherent fermion-hole excitations behave like bosons. This 
phenomenon is described by the Luttinger model. The main physical applications of this 
model occur in the context of the quantum Hall effect j3J |3] . However, the model is clearly 
also relevant to a trapped fermion system such as we have described. The field of fermion 
trapping is advancing very rapidly 13 El El , and it may be possible in the forseeable 
future to observe Luttinger-liquid-type behavior in trapped fermionic gasses. 

We consider a cylindrical trap with hard walls, of length L and axial radius b ^ L. 
The fermion modes in this trap are described by three quantum numbers — k, £, and rip- 
corresponding to the cylindrical coordinates (z, 6, p). k is the momentum along the axis of 
the trap (setting h — 1), while I is the angular momentum about this axis, and n p is the 
number of radial nodes in the wavefunction. The wavefunction is periodic the ^-direction 
with period L and vanishes at p = b. In the absence of other interactions, the energy of 
a trapped particle is 



where m is the fermion mass and is the n p -th root of the Bessel function Jg. For 
fixed k, the three lowest-lying levels have n p = and I = +1, 0, and —1. The energy gap 
for the two \£\ = 1 states is A = « Jff . 

Previously, we have considered the effects of virtual particles in radially and angularly 
excited states of the trap on a Luttinger liquid lying entirely in the ground state [TUJ. In 
this paper, we consider a different regime, in which the Fermi energy ep is larger than A, 
so that the lowest-lying i = ±1 states are filled with large numbers of particles, while all 
the higher excited states of the trap remain empty. We shall treat this system as consisting 
of three interacting Luttinger liquids and shall study the normal modes of excitation. For 
particular choices of couplings, there can exist extremely simple exact solutions of this 
model. 

In accordance with the Luttinger model |TTJ E] , we linearize the fermion spectrum 
and append an infinite negative energy Dirac sea to the physical system. If we fix the 
number of particles with each value of £, then we can take the noninteracting Hamiltonian 
(ignoring any additive constant) to be 
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where the p± operators are 
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The a (6 ) operators correspond to right- (left-) moving (spinless) fermions of angular 
momentum £ about the axis of the trap, and v F is the Fermi velocity for the Luttinger 
liquid of angular momentum £. (For quantities, such as the Fermi velocity, which depend 
only upon \t\, we shall never write a negative superscript; only "0" or "1" will be used.) 
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The p ± obey the Bose commutation 

relations \p e + (-q), P+(q')\ = [pt(q), p l L (-g')l = V^'lf and P-(l')\ = °> so the y 

create and annihilate phonons with energies \q\vp. 

We may also introduce fermion-fermion scattering terms which are bilinear in the p± 
operators. At low temperatures, nearly all the scattered particles lie near the Fermi sur- 
face, with momenta k F = mv F . We shall consider both forward scattering interactions 
of the types (k F ; kp) — > (kp; k F ) or (hp; —kp) — > (kp; —kp) and backward scattering 
(kp;—kp) — > (—k F ;k F ). The backward scattering is forbidden by momentum conserva- 
tion unless £ = ±£', and it is truly distinct from the forward scattering only if £ ^ £'. 
["Inelastic" scattering processes of the form (kp; —kp ) — > (k F ; —k F e ) with \£\ \£'\ are 
also allowed by angular and linear momentum conservation, but they lead to substan- 
tially different and more complicated effects, in part because they change the numbers of 
particles in the three Luttinger liquids. We shall ignore these interactions in the body of 
this paper, although a brief discussion of them is located in the Appendix.] 

We write the forward scattering terms as generalizations of the interactions in the 
Tomonaga-Luttinger model [131: 
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All the interaction strengths gu\\e\ are presumed to have the same underlying origin; how- 
ever, their values depend upon the radial density profiles of the two interacting particles, 
so they depend upon \£\ and \£'\. The g\i\\i>\ may most generally be functions of the 
momentum transfer q, but their large-g behavior is irrelevant to the low-energy physics, 
and we shall take the interactions to be momentum- independent. The symmetries of the 
matrix in (j3j) are physically very reasonable. The £ = states interact symmetrically 
with the £ = ±1 states, and the forward scattering interactions among \£\ = 1 states are 
independent of the sign of I. 
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We shall introduce the remaining Hamiltonian for backward scattering later, after we 
look a bit further at the forward-scattering effects. We may simplify the interactions by 
introducing symmetric and antisymmetric combinations of the \t\ — 1 operators: p± — 

73 (Pi* + P± X ) an d P± = 775 (Pi* — Pi 1 ) • This change of basis mixes states which are 
degenerate under H , and it simplifies H is to 
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The antisymmetric modes are completely decoupled from the others. 

The separation off of the antisymmetric \i\ = 1 modes is extremely fortuitous, because 
the same separation also occurs when backward scattering is included. The backscattering 
Hamiltonian is 

-"bs = y 2.^ fi'bsOfc Op a p+2k F +q°p-2k F -q- {') 

k,p,ql,l'=±l 

We have neglected the backscattering for the I = modes, because it may be rewritten as 
an effective forward scattering (i.e. as a contribution to goo), i^bs adds additional terms 
to the linear coupling matrix and introduces a sine-Gordon coupling jTJ HH HH] as well: 
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*> A is the boson field associated with the antisymmetric modes 
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The parameter a 1 is a short- distance cutoff, and v l F a l is the bandwidth of the I = ±1 
Luttinger liquids. (In [TUj, we found that a 1 = 1/kp for a system of trapped fermions 
such as we are considering.) 

Much of the remaining work we shall do on this model involves Bogoliubov transform- 
ing the operators in order to simplify the interactions, replacing the operators p± (q) with 
new operators p±{q) that obey the same commutation relations [TSj. Our transformations 
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are formally unitary only if the interactions decrease sufficiently quickly as functions of 
|g|. However, we have previously stated that the large-momentum behavior of the cou- 
plings is irrelevant to the low-energy physics, and for momenta above 1/c/, the model is 
unphysical in any case. So we shall continue to treat the couplings as constants. 

Performing an appropriate Bogoliubov transformation on the antisymmetric modes 
transforms the corresponding portion of the Hamiltonian into 
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where the p A and (f> A are the transformed operators and the corresponding field, and ua = 

J (vp) 2 — ((7bs/27r) 2 . For repulsive interactions (gbs > 0), the cosine term is renormalized 
to zero at long wavelengths. For < 0, the flow is not toward g^s = 0, but for a 
particular choice of couplings, the Hamiltonian is exactly solvable |17j . 

The self-coupled sine-Gordon field 4> A describes one set of fundamental modes of the 
system. We now turn our attention to the remaining modes. With another Bogoliubov 
transformation, the remainder of the Hamiltonian becomes 
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The transformed coupling is 
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We shall denote the matrix appearing in (jllj) by M . 

If the interactions g 00 , gn, and g hs conspire to produce a degeneracy, with u = u s = 
u, then M may be diagonalized easily. We again take symmetric and antisymmetric 
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form 
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One further Bogliubov transformtion gives the decoupled normal modes, with frequencies 
ui = \q\\/u 2 ± 2ug. 

In the presence of a m = u$ degeneracy, we find exact energy shifts containing terms 
linear in g. In the absence of this special degeneracy, we expect the lowest-order correction 
to the frequency to be 0(g 2 ). We may find this leading term by treating g perturbatively. 
Through a completely straightforward calculation, we find shifted frequencies 

1-2 

~ \q\uQ + \q\g — 2 ( 16 ) 

ii 11-2 2m /1r7 » 

u s ~ \q\us-\q\g— j ( 17 ) 

n - u s 

for the p° and p s modes, respectively. 

If we wish to write H" (for general Uq, us, and g) as a sum of noninteracting harmonic 
oscillator Hamiltonians, we must find a {7(2, 2) matrix that diagonalizes M. That is, we 
need to find a matrix V such that V'MV is diagonal, and V^I^^Y — h,2, where ^2,2 = 
diag [1 —1 1 —1]. [Our subsidiary condition on V — which makes it a member of U(2, 2), by 
definition — ensures that the transformation preserves the correct commutation relations. 
The Bogliubov transformation that led to the form (|TT|) for H" corresponds to an element 
of the block diagonal U(l, 1) x {7(1, 1) subgroup of {7(2, 2).] 

The {7(2, 2) diagonalization problem for M is equivalent to the ordinary diagonaliza- 
tion of 
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which is related to M by a phase rotation of the p_ terms. The positive eigenvalues of 
M' are exactly the normal mode frequencies, and the eigenvectors of M' give (upon an 
appropriate undoing of the phase rotation) the weights of the various p in the normal 
mode creation and annihilation operators. 

The positive eigenvalues of M' give normal mode frequencies of 



w = ^7|V M o + u s ± V( M o - u sY + I6g 2 u u s . (19) 

When Uq = us, we recover our earlier result, and for uq ^ us, the exact expression reduces 
to the perturbative result up to 0(g 2 ). The corresponding eigenvectors of M' are fairly 
complicated, but they may be found by elementary means. 

This gives a complete description of the modes of the three interacting Luttinger 
liquids. Two set of modes — which are linear combinations of p° and p s — describe linear 
harmonic fluids, with the two different wave speeds given by (|19|). The remaining modes 
are coupled through a sine-Gordon interaction ^SJ ^] . If we fix the number of particles 
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having each value of £, then, under the approximations of the Luttinger model, the p e 
operators generate a complete set of states, and the modes we have found are the only 
low-energy excitations for the system. Because the interactions iff s and ifbs are fairly 
general and physically quite reasonable, the spectrum we have derived should be highly 
relevant to the study of trapped fermions in situations where Luttinger liquid behavior is 
observable. 
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Appendix: Hamiltonian for "inelastic" scattering 

In the body of this paper, we have considered only those scattering processes which pre- 
serve the individuals particles' I values. However, there exists another type of scattering 
interaction, which does not have this property. Interactions of the type (k e F ; —kp £ ) — > 
(k F ; —kp £ ) with \£\ \£'\ are also allowed by linear and angular momentum conservation. 
In this appendix, we shall examine the effects of such terms. We shall find that they lead 
to interactions that are similar to sine-Gordon interactions but more complicated. 

We begin with an interaction Hamiltonian that describes the scattering process in 
question: 

H * = 7 E E 9in4b-_t +q a Jb \ +q + h.C, (20) 
k,p,q£=±l 

where "h.c." denotes the Hermitian conjugate. We may simplify this by Fourier trans- 
forming the a and b operators and expanding the resultant fermion operators in terms of 
the boson field. For example, a^ 1 may be written as [T5J 1201 ED 1221 EH] 

■ +1 - 1 f dxe- ikx iPl l {x) 
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v / 27ra r L A) 

The operator U+ 1 carries the fermion number of 1 but does not affect the bosonic state 
of the system, while <p +1 is the boson field corresponding to the I = +1 fermion states, 
and 9 +1 is the integral of the momentum conjugate to <p +1 . In general, the operators </> 
and 6 l have the forms 

-itf = JL ]T -e-^l/V^V + (p) + pi (P)] + in Nt+ + NL x (22) 
L p^O P L 
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where and N_ are the numbers of excess right- and left-moving fermions of angular 
momentum £, respectively. 

If we express all the operators in H in in this (bosonized) fashion, we get 



fjh 



S / L ^idx 2 dx 3 rfa;4e ifc(:C2 - a:i) e-^ (x2 - 2:i) e i9(:C4 - a:2) e ip(:E3 - S4) 



(2tt)%W fe; „ / ,_ ±J . 

The sums over fc, p, and q give ^-functions according to J2k e%kx = L6(x), and these 
simplify the expression substantially. We find that 
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Since ^ and 9 are Hermitian, this is equivalent to 
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This a new nonlinear interaction. It is similar to a sine-Gordon interaction, but it contains 
substantial additional complexity. The effects of this interaction would presumably need 
to be determined perturbatively. 
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